Spin dynamics in the ordered spin ice Tb2Sn207 
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Georaetrical frustration is a central challenge in 
contemporary condensed matter physics, a crucible 
favourable to the emergence of novel physics. The py- 
rochlore magnets, with rare earth magnetic moments lo- 
calized at the vertices of corner-sharing tetrahedra, play 
a prominent role in this field, with a rich variety of ex- 
otic ground states ranging from the "spin ices" Ho2Ti207 
and Dy2Ti207 to the "spin liquid" and "ordered spin 
ice" ground states in Tb2Ti207 and Tb2Sn2 07. Inelas- 
tic neutron scattering provides valuable information for 
understanding the nature of these ground states, shed- 
ding light on the crystal electric field (CEF) level scheme 
and on the interactions between magnetic moments. We 
have performed such measurements with unprecedented 
neutron fiux and energy resolution, in the "ordered spin 
ice" Tb2Sn2 07. We argue that a new interaction, which 
involves the spin lattice coupling through a low temper- 
ature distortion of the trigonal crystal field, is necessary 
to account for the data. 

In spin ice materials such as the pyrochlores IIo2Ti207 
and Dy2Ti207, the orientations of the Ising-like rare 
earth magnetic moments obey local constraints, the so- 
called ice rules, analogous to those which govern the hy- 
drogen vector displacements in water ice, and yield the 
same ground state entropy [1]. The ice rules do not en- 
force long range order, and many short range ordered 
configurations have the same ground state energy. Start- 
ing from these Ising configurations, a spin fiip which vi- 
olates the ice rules induces magnetic excitations akin to 
magnetic monopoles, propagating along Dirac strings [2 . 
The energy necessary to create a monopole is directly re- 
lated to the Ising character, namely to the crystal electric 
field (CEF) scheme of the rare earth ion considered. In 
the pyrochlore structure, the crystal field of Ho or Dy 
ions yields a ground doublet with strong Ising character 
separated from the first excited doublet by a large en- 
ergy gap A of 200 and 300 K respectively. Weakening 
the Ising character of the wave- functions and decreasing 
this energy gap, as it occurs in Tb2Ti207 and Tb2Sn207, 
where A 10-20 K, provides an extra degree of freedom 
in the spin ice structure by reducing the anisotropy en- 
ergy. Whereas it is not clear whether monopoles can still 
be defined in such a case, this degree of freedom allows 
new kinds of magnetic orders, possibly induced by ad- 
ditional perturbations, to be stabilized. Presently, the 
spin liquid state in Tb2Ti207 [3 and the ordered spin 
ice state in Tb2Sn207 [4 , already studied for some time. 



are still not fully understood [S] [6] . 

In this field, the ordered spin ice Tb2Sn2 07 offers a fas- 
cinating challenge. Recalling that an effective ferromag- 
netic first neighbour interaction is needed to derive the 
ice rules in the model spin ices, one can expect that, by 
decreasing the strength of the anisotropy, the short range 
ordered canonical spin ice, with moments along their lo- 
cal <111> axes (see Figjl] a), would transform into a 
long range ordered Heisenberg (isotropic) ferromagnet, 
with moments along <001> (see Figjl] b). The interme- 
diate state for finite anisotropy would be akin to a canted 
non-collinear ferromagnet. The validity of such a picture 
was checked by Monte-Carlo simulations [7 , where first 
neighbor ferromagnetic interactions in the pyrochlore lat- 
tice were combined with a uniaxial anisotropy term DS^ 
of variable magnitude. In this "soft spin ice" model, the 
evolution of the magnetic structure as the anisotropy is 
decreased can be monitored by the variation of the mo- 
ment angles with respect to their local anisotropy axes, or 
equivalently by the evolution of the reduced magnetisa- 
tion r = Mo/Msat along <001>. The ratio r varies from 
l/y/?> = 0.577 to 1 when going from Ising to Heisenberg 
behaviour. 

At first sight, Tb2Sn207 shows most of the charac- 
teristics of this simple model. The magnetic order pa- 
rameter deduced from the magnetic Bragg peak inten- 
sities [4 undergoes an S-shape rise as temperature de- 
creases, with a progressive onset below 1.3 K, an abrupt 
increase at Tn=0.87K, then a quasi-plateau below 0.6 K. 
The magnetic ground state reminds that of a canted fer- 
romagnet: the four tetrahedra of the pyrochlore cubic 
unit cell obey the ice rules in a coherent manner, yield- 
ing a non-zero global magnetization. For this reason, 
the magnetic structure found in Tb2Sn207 has been la- 
belled an "ordered spin ice" (OSI). However, a detailed 
inspection reveals features which do not match with the 
simple "soft spin ice" model described above [7 : i) the 
negative exchange contribution to the Curie- Weiss tem- 
perature [8 (6>cw — — 6K), together with the first mea- 
surements of the dispersion curves of the spin excitations 
[Hi [To], strongly suggests the presence of antiferromag- 
netic first neighbor interactions, which is difficult to rec- 
oncile with a ferromagnetic coupling; ii) the spin struc- 
ture deduced from a symmetry analysis of the neutron 
diffraction data and the derived canting angle of the Tb 
moments show that they are tilted away from the local 
<111> anisotropy axis, in the direction of the <110> 
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(a) : r = l/VS (c) : r < l/VS 




FIG. 1: (color on line) Different magnetic configurations on 
a tetrahedron : (a) "Two-in two out" configuration typical of 
the canonical "spin ice" ground state with magnetic moments 
along the local <111> directions; (b) Standard ordered fer- 
romagnetic configuration; (c) "ordered spin ice" ground state 
as found in Tb2Sn207. 



axis, which leads to a spontaneous reduced magnetiza- 
tion r 0.27 (see Figjljc) well below the range predicted 
by the soft spin ice model of RefEl 

To study this question, we performed detailed mea- 
surements of the dynamical structure factor S{Q^uo) in 
Tb2Sn207 by means of inelastic neutron scattering ex- 
periments (see Methods). We used the time of flight 
spectrometer IN5, installed on a cold neutron source at 
the Institut Laue Langevin (ILL), in its recent high in- 
tensity and high resolution set-up Our precise data 
are in agreement with previous measurements [10 ; they 
yield the possibility to test a fine scheme of interactions, 
taking both the energy and intensity dependence of the 
magnetic scattering into account. Our measurements un- 
veil strong similarities with the titanate Tb2Ti2 07 which, 
despite a similar CEF scheme, remains in a spin liquid 
state with no magnetic order down to T=0 [3]. Based 
on experimental evidences, we recently proposed that a 
strong 4f electron - strain interaction [T2J, resulting in a 
small "hidden" distortion developing at low temperature 
which breaks the local trigonal symmetry, may be a new 
promising route to understand this puzzling material ^ . 
The distortion induces a two singlet ground state, and we 
found that the spin liquid behavior arises from quantum 
spin fluctuations between the singlets. In the present case 
of Tb2Sn207, we assume that the same physics is also 
relevant and examine the influence of a similar "hidden" 
distortion. Indeed, to account correctly for the neutron 
data, and especially for the strong intensity of the lowest 
energy excitation, we show that one needs to consider 
not only the exact trigonal crystal field scheme [8 and 
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FIG. 2: (color on line) : Mapping of the powder average 
S{Q,uj) in Tb2Sn207 measured at 1.3 K (top) and 65 mK 
(bottom). These data were taken on the IN5 spectrometer 
installed at ILL, France, with an incident wavelength A= 5 A. 



the interionic interaction energy term, but also the dis- 
tortion from trigonal symmetry. Our model reproduces 
both the dynamic structure factor S{Q,uj) and the static 
properties in the magnetic phase, i.e. the ordered spin 
ice structure with the correct canting angle and modulus 
of the spontaneous Tb moment. 

The upper panel of Fig j2] shows a mapping of the pow- 
der average S{Q^uj) taken at T= 1.3 K in the paramag- 
netic regime, just above the transition. One observes a 
quasielastic signal at energies less than 0.4 meV, a gap 
in the 0.4-0.9 meV region and a broad energy band in 
the 0.9-1.6 meV region. This band is strongly modulated 
with a minimum around 1.1 A~^. All these features re- 
call those observed in the spin liquid Tb2Ti207 [8l fT3lfT4] . 
The modulation is due to antiferromagnetic correlations 
of Tb^+ moments over first neighbour distances within a 
given tetrahedron, so that the energy of the modulation 
reaches a minimum nearby the liquid-like peak, whose 
position is related to the first neighbor distance. 

At 65 mK (lower panel of Figj2|, the large quasielastic 
band splits into two parts. The band centered at zero 
energy becomes limited to a small region below 0.1 meV, 
where the traces of the magnetic and nuclear Bragg peaks 
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FIG. 3: (color on line): (a) Energy scans for a wavevector 
transfer Q=0.5 and several temperatures in the OSI 

phase. These cuts are obtained from the map by integrat- 
ing over a ±0.04 wide band around Q. (b) Temperature 
evolution of the peak position around 0.4 me V along with its 
energy width, deduced from phenomeno logical fits performed 
at several temperatures. 



can be detected. Well separated from it, a low energy in- 
elastic line, which is weakly dispersive, is clearly observed 
around 0.4 meV. This is further illustrated in Figjsja), 
which shows several scans as a function of energy trans- 
fer and fixed Q in the temperature range 65mK-l.lK. 
At the base temperature and up to 0.7 K, the low en- 
ergy inelastic line is well resolved. The high energy band 
also undergoes some changes and presents more structure 
than in the paramagnetic phase. Figjsjb) shows the tem- 
perature evolution of the peak position around 0.4 meV 
along with its energy width. While the peak is clearly 
observed in the ordered phase, it becomes overdamped 
in the paramagnetic regime. Understanding the origin of 
this low energy transition appearing in the OSI phase is 
the main goal of the present study, and it will in turn 
lead to a realistic understanding of the OSI phase. 

To analyze these experiments, one must go beyond 
the "soft spin ice" model. First, one must consider the 
actual CEF hamiltonian acting on a Tb^+ ion (J=6, 
gj=i/2): ^CEF = E„,„^ B^O^, where the are the 



usual Stevens operators equivalent. Second, besides the 
isotropic short range exchange J7 of a Tb ion with its z=6 
nearest neighbours, the infinite range dipolar interaction 
must be included (see Methods) . In the nominal trigonal 
symmetry, the parameters and the precise nature of 
the Tb^+ CEF wave functions and energies have 
been investigated in detail in Ref.8, The Tb^+ ground 
state is a degenerate doublet, close to: 

|^±)~|J = 6;J, = ±5), (1) 

while the first excited state is another doublet — 
I J = 6; J;, = ±4) with energy Ee - Eg = A 1.4meVc^ 
16 K. The dynamic structure factor S{Q^uj) observed in 
the paramagnetic phase is thus readily interpreted: the 
width of the quasielastic band reflects fluctuations of 
the Tb moments inside the ground state doublet IV^^), 
whereas the broad band at finite energy corresponds to 
spin excitations from the GS doublet to the first excited 
doublet modulated by exchange. 

The main evolution of the spectra between 1.3 K and 
65 mK, i.e. the appearance of the inelastic line at 
0.4 meV, could at first sight be explained by the triv- 
ial effect of the molecular field experienced by a Tb ion 
in the OSI phase. At the base temperature, when only 
the lowest energy level is populated, this 0.4 meV tran- 
sition would result from the exchange/dipole splitting of 
the GS doublet, whereas the changes in the high energy 
band would arise from that of the first excited doublet. 
The new GS eigen-functions in the OSI phase, describ- 
ing the exchange split doublet, do not strongly depart 
from the zero field states ([T]) because of the modest mag- 
nitude of the molecular field h (0.5-1 T, see below). In- 
deed, admixtures of the excited doublet IV^^), which oc- 
cur through the transverse component h± (J+ and J_ 
operators), are of order h±/A and remain small. This 
point is of fundamental importance for our purpose be- 
cause, for non-Kramers ions like Tb^+, the operator J, 
which is odd with respect to time reversal, has no ma- 
trix elements between the time conjugate states of the 
doublet 

{^t\J\^;) = {^t\J\^;) = o, (2) 

and this also holds to a very good approximation for the 
exchange split states. Since the spectral weight of a CEF 
transition between states |?/^^) and li/j^) observed by neu- 
tron scattering is proportional to \{iIj^\J\iIju)\'^ , the tran- 
sition between the molecular field split states of the GS 
doublet is expected to have a very small intensity. Hence, 
the corresponding line should be extinct, contrary to ex- 
perimental observation. 

To put these arguments on a more quantitative ground, 
we numerically diagonalise the total hamiltonian in the 
mean field approximation and compute S{Q^uj) in the 
standard RPA approximation (see Methods), giving ac- 
cess to the dispersion of the collective spin excitations. 
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FIG. 4: (color on line) Calculated powder average <S(Q,cj) 
at T = 65 mK in the RPA approximation for J' = — 0.09K, 
using the trigonal crystal field of Tb2Sn207. 



In Tb2Sn207, an antiferromagnetic exchange integral 
J' = — 0.09K can be inferred from the analysis of the 
Curie- Weiss law (see Methods). Figj4] shows that, at 
65 mK, 5(Q,cj) consists in a set of crystal field tran- 
sitions around 1.5 me V slightly modulated by the ex- 
change/dipole coupling between magnetic moments. In- 
terestingly, despite the ferromagnetic-like nature of the 
ground state, the minimum of the calculated dispersion 
curve does not occur at the ferromagnetic zone center 
(Q=0), but at its boundary, i.e. it is in line with the an- 
tiferromagnetic exchange integral. However, as expected, 
a comparison between Fig j4] and the lower panel of Fig|2] 
shows that the model is not able to capture the intensity 
of the low energy transition observed at 0.4 meV. 

Explaining the data thus demands the introduction of 
a new interaction term that would imply a major change 
in the CEF wave functions. At this step, it is worth 
emphasizing that, for a non-Kramers ion, a degenerate 
doublet is very unstable against any crystal field per- 
turbation. Several arguments suggest that a tetragonal 
distortion is at play in the sibling material Tb2Ti207: 
there exists clues to the presence of a T Jahn- Teller 
transition \TT and of cubic-to-tetragonal structural fiuc- 
tuations below 20 K [15 . The distortion is likely to be 
at the origin of the spin liquid state of this material [6]. 
We assume here that such a distortion is also present 
in Tb2Sn207 and we examine the consequences of this 
hypothesis on the physical properties. With its principal 
axis along the cubic [001] axis, the distortion is described 
by the following crystal field operator equivalent, in the 
local frame with the <111> axis as z-axis 116 : 



Let's examine first the changes introduced by the distor- 
tion on the ground state in the paramagnetic phase. As 
in the case of Tb2Ti207, the principal effect is a degener- 
acy lifting of the ground doublet (and of the first excited 
doublet) and a deep modification of its wave-functions. 
Whatever the nature of the distortion, it can easily be 
shown that the new eigenstates become entangled sin- 
glet states of symmetric and antisymmetric type which, 
for the ground state, are close to: 



V2 



]J = 6;J, = 5)±|J = 6;J, 



-5)]. (4) 
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[2Jl -^Jl^V2 {J^J, + J, J,)]. (3) 



These states are such that {iIjs\J\^I^s) = ('0o|-^|'0o)=O, but, 
by contrast with the non-perturbed states ([T]), the cross 
matrix element (il^alJzli^s) is non- vanishing and large. 
For distortion amplitudes Dq similar to those derived 
in Tb2Ti2 07 (0.25 K [6 ), the energy separation between 
the singlets amounts to a few 0.1 meV (a few K). There- 
fore, we expect a new low energy transition in the CEF 
scheme from l?/^^) to with a large spectral weight. 

In the OSI phase, for the same reason as mentioned 
above, the eigen-functions do not depart strongly from 
the singlets Q, resulting in a sizeable intensity low en- 
ergy transition, as observed experimentally. The inten- 
sity map calculated in the presence of the distortion for 
J7 = — 0.09K is represented in Fig|5j in the param- 
agnetic phase (upper panel) and in the ordered phase 
(lower panel), to be compared with Figj2] Fig J6] shows 
ID energy cuts (experimental and calculated) for two 
Q values at 65 mK. As expected, at the base tempera- 
ture, a sizeable intensity for the low energy line is re- 
covered, in agreement with experiment, and its posi- 
tion (0.4 meV) is correctly reproduced with a distortion 
strength Dq = 0.2 K. Its intensity, however, is lower than 
observed. Preliminary calculations show that the agree- 
ment with experiment can be improved by introducing 
anisotropic exchange, but this will not be discussed here. 
The Q-dependence of the structure in the upper broad 
band is also rather well captured by the calculation, at a 
semi-quantitative level reproducing the inversion of the 
energy asymmetry between (5=0.5 and Q=1.0A~^ . In 
the paramagnetic phase, for Dq = 0.2 K, the splitting 
between the two singlets should be J 0.17 meV. The 
corresponding inelastic line is thus expected to merge 
into the quasi-elastic peak, which extends up to 0.4 meV 
at 1.3 K, and this explains why it cannot be clearly ob- 
served in the experiment. The above S value is compati- 
ble with the thermal evolution of the energy of the inelas- 
tic line shown in the lower panel of Figjsj which reaches 
0.17 me V at 1.2 K and becomes overdamped above this 
temperature. 

It seems therefore clear that the introduction of a 
tetragonal distortion greatly improves the agreement be- 
tween the calculated dynamic structure factor and the 
inelastic neutron data. We turn eventually to the static 
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FIG. 6: (color on line) Experimental (left) and calculated 
(right) inelastic spectra at 65 mK, for Q — 0.5 and 1.5A~^, 
in Tb2Sn2 07. The simulations were performed with J — 
— 0.09K and Dq=0.2K. A constant background was added 
to the calculated profiles, and their overall intensity was ar- 
bitrarily scaled to reproduce the intensity of the upper band 
at Q = lA"^ 
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FIG. 5: (color on line) Calculated powder average S{Q^uj) 
with a tetragonal distortion of the trigonal CEF {Dq = 
0.2 K), for J — — 0.09K, in the paramagnetic phase (T = 
2.5 K, upper panel) and in the ordered spin ice phase (T = 
65 mK, lower panel). 



properties of the OSI phase in Tb2Sn207 by comput- 
ing the magnetic phase diagram (j7,Tc) using our 4-site 
self-consistent model (see Methods). Without distortion, 
as expected, large negative (AF) J values {J < Jc = 
—0.22 K) stabilize a non-collinear Neel ground state, the 
so-called "all in - all out" antiferromagnetic structure. 
Above this critical value, the ground state is the ordered 
spin ice phase (OSI) with the same characteristics as 
observed by neutron diffraction in Tb2Sn207 [4 . The 
existence of this ferromagnetic-like phase with an AF 
exchange integral is fully consistent with the now well 
established idea that, in the presence of dipolar interac- 
tions, a spin ice ground state can be stabilized for moder- 
ate AF exchange [17 . This result is also consistent with 
the "soft dipolar spin ice" model [TS] recently proposed 
to account for the ground state of Tb2Sn207. Our mean 
field model reproduces a peculiarity of the OSI phase, al- 
ready mentioned in the introductive section, that puts it 
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FIG. 7: (color on line) Ordered spin ice part of the {J,Tc) 
phase diagram computed using the crystal field interaction in 
Tb2Sn2 07 in the presence of the tetragonal distortion (|3]), for 
different values of the Dq parameter. 



outside the reach of the "soft spin ice" model of Ref,7: the 
magnetic moment is found to tilt away from the <111> 
axis, in the direction of the <110> axis, or equivalently 
the reduced magnetisation r is lower than l/\/3. Figure 
[T] focuses on the OSI part of the phase diagram, com- 
puted for different amplitudes Dq of the distortion. For 
Dq < 0.25 K, the critical temperature Tq shows a regular 
variation with J'. For the parameter values characteristic 
of Tb2Sn207: J = -0.09 K and Dq 0.2 K, we find a 
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transition temperature Tq — 2.35K, higher than the ex- 
perimental value 0.9 K [4J. We think this difference is 
due to a possible first order character of the transition 
at 0.9 K, suggested by the abrupt S-shape variation of 
the magnetic Bragg intensity near Tq. The saturated 
Tb moment is 6.15 /iB, close to the experimental value 
5.9(1) /iB, and the angle between the moment and the 
<111> axis is 15°, also close to the experimental value 
13°. The total molecular field (exchange + dipolar) is 
h = 0.78 T and it makes an angle of 55° with the local 
<111> axis. 

For higher distortion amplitudes, a dip appears in the 
Tc{J) curve for low \J\ values, and above a distortion 
amplitude Dq 0.4 K, one observes a strong depletion 
of the ordered spin ice critical temperature, which goes 
to zero in a narrow range around J' :^ — 0.03K. This 
phenomenon is reminiscent of the problem of the sin- 
glet ground state in non-Kramers rare earth systems [19] , 
where no magnetic order occurs down to T = if the ex- 
change integral is below a critical value proportional to 
the gap separating the singlet from the excited states. In 
a recent work [6 , we argue that this case is realized in 
the spin liquid Tb2Ti207. Likewise, a sizeable intensity 
low energy transition appears in the neutron scattering 
spectra at low temperature, both in zero field and in the 
presence of an external field, which cannot be indexed in 
the trigonal CEF/Zeeman level scheme [6l[T6]. Our pic- 
ture points to the essential role of the precise CEF scheme 
of Tb2Sn207 and Tb2Ti207, as well as to the relevance 
of a new interaction, namely the 4f electron - strain cou- 
pling, to understand the properties of these compounds. 
In this sense, our model is very different from the many 
body approach described in [5], which emphasizes the 
importance of virtual crystal field excitations. 

In summary, we have measured low energy spin ex- 
citation spectra in the "ordered spin ice" pyrochlore 
Tb2Sn207 with a very high resolution. In the frame 
of a mean field calculation, we have shown that the 
ferromagnetic-like ordered spin ice structure is the 
ground state in this material for moderate antiferromag- 
netic exchange and taking into account dipolar interac- 
tions between Tb moments. We also find that the intro- 
duction of a small lattice distortion from trigonal sym- 
metry is necessary to explain most features of the low 
temperature dynamical structure factor, in particular the 
existence of a low energy transition in the ordered spin ice 
phase. Indeed, due to a specific property of non-Kramers 
doublets, the neutron scattering intensity between trig- 
onal symmetry states split by exchange is vanishingly 
small. The perturbation brought up by the distortion 
results in "entangled" tunnel-like wave- functions for the 
two lowest levels, yielding a sizeable intensity inelastic 
line, as observed. Our model reproduces the details of 
the magnetic structure in Tb2Sn2 07, and our results 
strongly suggest that the 4f electron - strain coupling 
is a key feature to analyze the quantum spin ice states in 



terbium pyrochlores. The physics underlying the "oppo- 
site" behaviours of the sibling materials Tb2Ti207 and 
Tb2Sn207, i.e. respectively the spin liquid and the or- 
dered spin ice phase as ground state, appears to be de- 
scribed by the same model, but with slightly different 
parameters. 
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METHODS 
Experimental details 

A powder sample of 2.2 g, the same as used in ReflU 
was inserted in a rectangular aluminium sample holder to 
minimize the absorption corrections, then placed into a 
dilution insert inside a helium cryostat. The temperature 
was varied from 65 mK to 1.6 K. For the measurements on 
the IN5 spectrometer, we used three incident wavelengths 
of neutron beam: A = 5, 8 and 10 A, yielding respectively 
an energy resolution of 90, 20 and 15 /ieV, at full width 
half maximum (FWHM) at the elastic peak. Spectra 
were corrected from the background of the empty sam- 
ple holder and from neutron absorption. The detector 
efficiency was taken into account by measuring a vana- 
dium plate of the same dimensions as the sample. The 
energy resolution was determined with a vanadium foil 
inserted in a cylindrical sample holder. 

Mean field calculation 

The total mean field Hamiltonian for a Tb^+ ion at a 
given site i reads: 

z oo 

U = UcEF + - MJ2 ^'.i + E ^diPi,,)(^) (5) 

J = l 3 

where Jij = J denotes the isotropic short range ex- 
change of a Tb moment with its z=6 nearest neighbours, 
and Jdip corresponds to the infinite range dipolar inter- 
action calculated using Ewald summation In defin- 
ing the neighborhood of a given moment, it is worth 
emphasizing that we assume a ^=0 magnetic structure 
in the FdSm space group with face centered cubic (fee) 
symmetry. In other words, the 4 Tb moments of a given 
tetrahedron may be different, but the spin configurations 
on tetrahedra connected by fee lattice translations are 
the same. We consider the 4 moments at the vertices of a 
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tetrahedron, diagonalize 1-L in the local frame for each site 
to determine the energies and wave functions IV^^), 
and perform a self-consistent calculation to obtain the 
Tb^+ magnetic moments —gjfiB{Ji) at a given temper- 
ature. The only parameters are the exchange integral J7, 
the lattice parameter a and the distortion magnitude Dq . 
The AF exchange integral in TbsSnsOy is J = -0.09 K, 
a value derived from the molecular field constant A ob- 
tained in ReflHl using the relationship: 



mean field approximation. 



denote the occupation 



J 



6 kB 



(6) 



RPA calculations 



In the RPA approximation, S{Q,uj) is given by the 
spin-spin susceptibility : 

m,n 

2^ Ua/3 j Xm,ni^), 

a,l3=x,y,z 

where fm denote the position of spin at site m. The 3x3 
tensor Xij is solution of the equation: 



The tensor is the 3x3 bare susceptibility, whose ele- 
ments ^^"^ are given by : 



(3 



of level /i, = (exp — £^^//cbT)/Z. The Wi^i, are 
vectors iW^^^^W^^^^^Wl^^) corresponding to the ma- 
trix elements of the total angular momentum Ji = 
Y.n V '^i^v\^i^)hl^iv\ in the eigen-basis. 



(7) 
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